Statistics of the contact network in frictional and frictionless granular packings 
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04 ■ Simulated granular packings with different particle friction coefficient fi are examined. The dis- 

tribution of the particle-particle and particle-wall normal and tangential contact forces P(f) are 
computed and compared with existing experimental data. Here / = F/F is the contact force F nor- 
malized by the average value F. P(f) exhibits exponential-like decay at large forces, a plateau/peak 
near f — 1, with additional features at forces smaller than the average that depend on fi. Com- 
putations of the force-force spatial distribution function and the contact point radial distribution 
function indicate that correlations between forces are only weakly dependent on friction and decay 
rapidly beyond approximately three particle diameters. Distributions of the particle-particle contact 
angles show that the contact network is not isotropic and only weakly dependent on friction. High 
force-bearing structures, or force chains, do not play a dominant role in these three dimensional, 
unloaded packings. 
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I. INTRODUCTION 



The properties of granular materials, even static sandpiles, continue to perplex engineers and physicists alike [jp. 
Processing of granular materials play a central role in the pharmaceuticals industry as well as engineering communities, 
such as ceramic component design. For example, one may wish to evenly distribute the ingredients in a tablet or 
pill, or reduce the likelihood of component failure. However as yet there is no clear indication of how the individual 
particle properties determine the final state of the system. 

Although packings of frictionless, monodisperse, cohesionless, hard spheres have been well studied little is 
known about the effect of including particle friction. Recent discrete clement simulations of granular materials, where 
packings were generated for particles for different static coefficients of friction (i between individual particle pairs, 
showed that the local particle coordination of the packing varied strongly as a function of friction || . From a different 
perspective, experimental studies of static granular assemblies have shown many interesting facets of the stress state 
. of these systems. One method of analysis appears to dominate in describing the statistics of granular packings: 
\Q ' computations of the distributions of normal forces are de rigeur. 

Experimental studies on granular packings use a variety of techniques to measure the distribution of normal contact 
forces P(f) between particles and container walls where / = F/F - all measured normal forces F are normalized 

with respect to the average force F. The Chicago group 0-|7j utilized carbon paper to measure P(f) at the base and 
J ■ sides of a cylindrical container packed with glass spheres with a normal load applied at the top of the packing. Forces 
several times the average force were observed, with resolution down to the weight of a few particles. Blair et al. || 
measured P(f) for amorphous and ordered granular packings for particles with different values of /i, which varied by 
a factor of approximately three. These experiments demonstrated that P(f) is indiscriminate towards the effects of 
particle friction and structure of the packing, and the general form of P(f) remained robust within the resolution of 
'■^ ' the experiment. 

In a different experimental set up, L</>voll et al. Q used a pressure transducer device to measure P(f) at the bottom 
of an unloaded granular packing under its own weight, on a fixed substrate of particles glued to the supporting base. 
This experiment was able to resolve forces down to the weight of a few grain masses and showed that the spatial 
distribution of contact forces were correlated over a few particle diameters. Using a novel modification of the carbon 
paper technique, Tsoungui et al. || actually measured P(f) inside the bulk of a 2D packing. Despite the poorer 
statistics of this study, the results agreed well with Blair et al. M and Lc^voll et al. m. Experimental studies on static 
granular packings show that P(f) exhibits several generic features; an approximately exponential tail at large / and 
a plateau or peak near / ~ 1. 

Computational studies, such as contact dynamics or molecular dynamics, of compressed packings provide good 
comparison with the experimental data |I^,|Tl|. However, as yet, there has been no systematic study of the effects 
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of particle friction on the force distributions within a granular assembly. Here we show how the effects of friction 
change the behaviour of P(f) in the small force region but only weakly affect the large-/ region. We show that the 
local contact geometry of the packing is not isotropic and only weakly influenced by friction. We also discuss aspects 
of the force network whereby high force-bearing structures, or force chains, do not seem to be a dominant feature of 
these unloaded packings. 

We computed P(f) in the bulk of various packings (which is presently inaccessible in 3D experiments) that had 
settled onto either a rough bed or a planar base. We compared these results with P(f) for particles in contact with the 
flat base (similar to experiment) of a periodic packing and with the P(f) generated at the side walls of a cylindrical 
packing. We resolve the components of the contact force that are normal (n) and tangential (t) to the line of centres 
between two particles in contact. 

In the next section we briefly describe the model, though a more thorough description of the technique is available 
elsewhere p||l^. In section [II, we present results for the force distributions, force correlations, and the contact 
geometry. We also discuss some aspects of the force network with respect to a force cut-off scheme, highlighting some 
pros and cons of this method. In section IV we summarize and conclude this work. 



II. MODEL AND METHOD 



We performed three dimensional (3D) molecular dynamics simulations with N monodisperse, cohesionless, inelastic 
spheres that interact only on contact via a Hooke (linear) spring or a Hertz contact law and static friction Jl3| , p^| . 
Contacting particles i and j positioned at and Vj experience a relative normal compression S = |ry — d\, where 
r ij = r i — r j ■ which results in a force = F„ + F t . The normal and tangential contact forces are given by, 

(7TI \ 
knSllij - y7„V„J , (1) 



F t = f(6/d) (-fc t As t - | 7t v t ) , (2) 

where = r^/r^j, with r,j = |r^- 1 , v„ and v t are the normal and tangential components of the relative surface 
velocity, and fc n ,t and 7 Wj t are elastic and viscoelastic constants respectively. f(x) — 1 for Hookean springs and 
f(x) — yfx for Hertzian contacts. As t is the elastic tangential displacement between spheres, obtained by integrating 
surface relative velocities during elastic deformation of the contact. The magnitude of As t is truncated as necessary 
to satisfy a local Coulomb yield criterion F t < fiF n , where F t = |F t | and F n = \F n \, and /i is the particle-particle 
friction coefficient. For the present simulations we set k n = 2 • 10 5 mg/d, k t — jk n , 7„ = 50^/gJd. For Hookean 
springs we set 7t = while for Hertzian springs, 74 = 7n . For Hookean springs the coefficient of restitution e n> t, is 
related to 7 „ it through, 

e n ,t = exp(- 7lM i co //2), 

where the collision time t co i is determined by the contact frequency between two particles. For the parameters chosen, 
e n = 0.88 for Hookean springs. For Hertzian contacts e is velocity dependent Jlq] . We chose a time-step St = 10 _4 r, 
where r = y 1 dj ' g and g is the gravitational acceleration. 

Amorphous packings (with packing fraction cf> w 0.60) were generated by allowing an initially dilute system to settle 
under gravity acting in the vertical direction. Particles settled onto a bottom wall that was either a planar base or a 
bumpy bed of particles frozen into a close packed random configuration. This process was run until the kinetic energy 
of the system was much smaller than the potential energy || . The base had the same frictional and elastic properties 
as the particles being poured. 

Most of our results are for packings that are spatially periodic in the horizontal plane, i.e. we ignored the effects 
of sidewalls. Because of this, the pressure in a packing does not saturate with depth. Therefore, to make a direct 
comparison with experiment, our definition of the average force f n j = Fn,t/ F(z) n ,t, was normalized by F(z) n t , the 
average contact force at a depth z in the packing. The generation of these packings is fully discussed in Ref. ||. We 
also compared results for packings poured into a cylindrical container with 'flat' walls and the same properties as the 
particles. In this case, there is no need for depth-average normalisation, as the walls carry a significant fraction of the 
weight of the system . Results for the periodic packings with depth-average normalisation are consistent with the 
cylindrically confined packings and experiment. As a consequence, the depth-average normalisation proves to be the 
correct method for dealing with periodic packings. 
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III. RESULTS 



A. Force Distributions 

Force distributions in all granular packings exhibit several general features. Measurements of the distribution of 
normal contact forces P(f n ), for granular packings that are either free-standing under the influence of gravity j^] (as 
we simulate here) , confined packings that have been loaded (as in experiments) , or axially compressed systems 

(as in previous simulation studies and experiment) [[l0||l7|Jl8f| , as well as a lattice model ||, all purport exponential 
tails in P(f n ) at large forces (typically for f n > 1). Mueth et al. [Q used an empirical fit to their experimental data 
of the form, 

P(fn) = a(l-be-£)e- f, f«, (3) 

and found a — 3.0, b = 0.75, and (3 = 1.5 ± 0.1 for loaded glass spheres confined in a cylindrical container. 

In Fig. [l] we show our computations of the force distributions for the normal contact force /„ for different systems. 
In Fig. |l](a) we see that the form of P(f n ) is the same for both Hookean or Hertzian contact force laws. Varying 
the system size has no effect (other than improving the statistics of the data) . Similarly, in Fig. 0(b) computations 
of P(fn) m the bulk of a periodic or confined system, at the base of the periodic system, or at the sidewalk of the 
cylinder are indistinguishable. Recent 2D simulations have shown that P(f n ) at the base can depend on the properties 
and geometry of the base (T^|. Computations of P(f n ) for those particles in contact with the flat base and at the 
side walls also show the generic form seen in the other data although the statistics here are poor due to the number 
of contacts in the plane (sa 10 4 ) compared with the number of particle-particle contacts in the bulk (ss 10 5 — 10 6 ). 
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FIG. 1. Distribution of normal contact forces P(fn) for packings of N monodisperse spheres of diameter d and particle 
friction coefficient fj, — 0.5. (a) Comparison between a spatially periodic Hertzian packing with a square base of dimensions 
A = 20d x 20d and a Hookean packing with A = 40<f x 40d. System sizes are indicated in legend, (b) Comparison between 
two Hookean packings, one a spatially periodic system with N = 128000 and A = 20d x 20d, the other a confined, cylindrical 
packing of diameter D — 20d and N = 50000. 

To compare with existing experimental data, we fit Eq. ||to our data for the largest system. We show this comparison 
in Fig. ||. The P(f n ) computed over all contact forces is denoted by the solid circles in Fig. || with best fit parameters 
a = 2.55, b = 0.65, and f3 — 1.35, agreeing well with Eq. || up to / w 2, but falling off more quickly than Eq. || for 
large /. We find a better fit to Eq. || if we exclude the data for F n < mg, i.e. ignoring all data in the limit / — ► 0, 
essentially mimicking the finite resolution in experiment. This alters the average value such that our original data 
set has now been 'squeezed' together. We denote this data as the partial set in Fig. ^. The fit to Eq. || with a = 3.1, 
b = 0.78, and (3 = 1.55, is much better than when data for small forces is included. Our simulation data is in good 
qualitative and quantitative agreement with previous experimental results B and similar to Radjai et al. [jTof. 
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FIG. 2. Distribution of normal contact forces P(fn) for Hookean packings of N = 128000 monodisperse spheres and fi = 0.5, 
on a fiat base of dimensions 40d x 40d. The full P(fn) (solid circles) includes normal forces for all contacting particles and we 
fit to Eq. ^ (solid line) using a — 2.55, b = 0.65, and f3 = 1.32. For the partial P{fn) (open circles) we have excluded all forces 
less than the weight of one grain and renormalized /, finding a better fit to Eq. ^| with a = 3.1, b = 0.78, and f3 = 1.55. We 
have arbitrarily shifted the curve for the partial P(f n ) for clarity. 

The empirical fit of Eq. || is poor for the total bulk P(f n ) over a large range of the data and we only achieve 
agreement by renormalising our data, using the partial data set in Fig. || We also note that on closer inspection of 
existing simulation and experimental data, whether the tails of P(f n ) are truly exponential or not is questionable and 
may be an indication of the averaging technique used in computational studies [ pp| . 

The distribution of tangential forces P(ft) is shown in Fig. [| In comparison with the normal forces, P(ft) decays 
more slowly than P(/„). Fitting Eq. ||to the bulk data for the largest system (N — 128000), we find good agreement 
with a = 2.5, b = 0.7, and = 1.4. 
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FIG. 3. Distribution of tangential contact forces P(ft) for packings of N monodisperse spheres of diameter d, with particle 
friction coefficient /i = 0.5. System size is shown in legend. The line is fit to Eq. ^ for the largest system. 

While there is clearly some agreement on the behaviour of P(f n ) for large /„, the characteristic nature of the small 
force region of P(fn) remains in dispute. Experimental data show that P(f n ) approaches a finite value as f n — * 0. 
However, some numerical works have suggested that P(f n ) ~ * f° r small /„ [^j]. In Fig. [J we show the small force 
region of P(f n ) for packings identically prepared but with different particle friction coefficients. We do not show the 
full P(fn) curve as friction only weakly influences the behaviour of the large-/ region. However, our large system 
size data suggests that the exponential tail becomes slightly steeper with decreasing friction, i.e. (3 increases as /i 
decreases. The defining feature of these packings is that for purely frictionless systems, P(f n ) shows a well-defined 
peak at small forces, while for /i > 0, P{f n ) has & n upturn at very small forces. The amplitude of this upturn increases 
with increasing friction coefficient. 
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FIG. 4. P(fn) at small forces for packings with different particle friction coefficient /i. Frictionless packings (/i = 0) exhibit 
a well-defined peak in P(fn) near /„ = 1, whereas even for low frictional packings, an upturn appears in P(f n ) at very small 
forces. The amplitude of this upturn increases with increasing friction coefficient and the position of the peak also shifts to 
larger /„. Results are for Hookean packings with periodic boundary conditions in the horizontal plane, for N = 20000 on a 
rough, particle base with A = 20d x 20d. 



The Chicago group || studied different particle packings where /i varied by a factor of approximately three. 
Within the resolution of their experiment they did not find any systematic trend with friction. Because of the higher 
resolution in simulation, the following comments are relevant to such studies: the fraction of particle-particle contacts, 
or bonds, experiencing small forces increases with increasing fx even though the total number of contacts decreases 
with increasing /i || . Further study shows that the fraction of particles that are close to the Coulomb yield criterion 
F t ~ fiF n , i.e. those particle pairs that are most likely to undergo local plastic rearrangement, increase as \x — > 0. 
Indeed, we have previously reported H that frictionless packings are always isostatic, whereas frictional packings are 
hyperstatic and this may be related to the behaviour of P(f n ) at small /„. 

For completeness we show the corresponding distributions P(ft), for the tangential forces in Fig. [j| In this case, 
we do not find any significant systematic trend with \x. The role of /i in the determination of P{f n ,t), is subtle. In 
frictionless packings, P{f n ) does not show an upturn at small /„, therefore the generation of this upturn in frictional 
packings comes from the very presence of the frictional forces f t , which influence the nature of particle contacts such 
that P(fn) itself observes an upturn at small forces. 
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FIG. 5. P(ft) at small forces for packings with different particle friction coefficient ^t. Results are for a Hookean packing 
with periodic boundary conditions for N = 20000 and A = 20d x 20d. 



B. Force Correlations 



The spatial force-force correlation function J-(r) measures spatial correlations between forces separated by a distance 
r. We use the same definition as in Refs. IH, 

F( r ) = ^ 4- — 17\ — i — r~ ' ( 4 
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where is the distance between particle contacts i and j, and fi is the normalized contact force acting at contact 
i. In experiment, spatial force correlations can, at present, only be measured at container walls: the points of force 
contact coincide with particle contacts at the container surface lying in a 2D plane. The minimum separation between 
measurements in experiment is coincident with the particle size, r m j„ w d. In a 3D packing, contact forces transmitted 
at the points of particle-particle contacts are only restricted by excluded volume effects. For monodisperse spheres 
in 3D the minimum separation, r m j„ rs |. A locally four-particle pyramid configuration would give this minimum 
separation. Mueth et al. [Q found no evidence for spatial correlations between the contact forces within the resolution 
of their measurements. L^voll et al. using a different measuring technique, resolved their force data showing weak 
force correlations at the base of their packings which extend out to approximately five particle diameters || . This 
may only come about from the induced order of the sample at the container wall. 

Because of the restricted geometry of experimental measurements, we found it instructive to compare our compu- 
tations of T{r) for the normal contact forces within the bulk of amorphous packings, and sec how these might depend 
on (x. For comparison we also computed the correlation function between tangential contact forces for /i = 0.5. In 
Fig. ^ we present the spatial force correlation function for a frictionless packing (fi = 0) and a frictional packing 
(/i = 0.50). Within the bulk of the packing, forces are correlated, but only over short distances, extending to less than 
three particle diameters in the bulk, indicative of the diffuse nature of the force transmission network. However, the 
effect of friction on these correlations is very weak, with the frictional packing exhibiting only a very slight increase 
in local correlation. 
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FIG. 6. Spatial force-force correlation function T{r) for normal contact forces as defined in Eq. ^| computed within the bulk 
of a frictionless (circles) and frictional (/i = 0.5 - solid line) packing. The dotted line shows the corresponding correlation for 
the tangential forces when g, = 0.5. 

Similar to Mueth et al. {jj, in Fig. fj] we also show the radial distribution function g(r), between contact points inside 
the bulk of a frictionless {fi = 0) and a frictional packing (/i = 0.50). Clearly, the frictionless packing has a higher first 
peak, representative of the higher coordination of the frictionless packing compared with the frictional one and 
also local correlations between the positions of the contact points are stronger in the case of the zero friction packing 
indicating a more ordered distribution of contact points in the system. 
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FIG. 7. Radial distribution function g(r) of the contact points within the bulk of a frictionless (circles) and frictional (/j = 0.5 
- line) packing. 
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C. Contact Geometry 



We have so far shown that computations of P(f) for various particle parameters yield essentially the same data, 
except for small /. It is ironic then, that although the generic features of P(f) are a signature of the granularity of 
the system, it offers little distinctive information on the grain-level properties of the packing. Keeping in the spirit of 
particle pair information, in Fig. || we show the probability distributions for particle-particle contact angles defined 
in the local spherical coordinate system that bonds make with respect to the vertical (parallel to gravity direction). 
In Fig. U we compare packings with different /i (=0, 0.1, 0.5) and found that the distribution of contact angles has 
only a weak dependence on friction indicating that all the systems locally appear similar. In all cases, the majority 
of contact angles lie in the range 45° < 9 < 90°. 
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FIG. 8. Probability distribution functions P(9) for particle-particle contact angles, where 9 is denned in the local spherical 
coordinate system as the angle the particle pair makes with the vertical. 6 — is a vertical contact and 9 = 90° a horizontal 
contact. Packings with [i = 0, 0.1, 0.5 are shown. 

It is a simple exercise to further compute P(9) only between particles that carry a large force, i.e. to identify or 
distinguish between "weak" and "strong" forces, as Radjai and co-workers have done for compressed systems plp2| . 
In Fig. [| we compare P(9) computed between all contacting particle pairs and P{6) computed for the subset of 
particles in contact whose normal contact force /„ > fcut: where f cu t is some given threshold value. Here we set 
fcut = 2.0, i.e. all particles whose normal contact force is greater than the twice the average. Resolving the contact 
angle distribution according to a force cut-off as in Fig. |p reveals that high force-bearing clusters are more directional 
and the anisotropy grows with increasing f cut (not shown here). 
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FIG. 9. Distribution of contact angles P(9) of particles in contact for a packing with n = 0.5. We distinguish between P(9) 
computed between all pairs in contact (solid line), and a sub-set of particle pairs whose contact force is greater than some 
cut-off threshold /cut (dashed line). Here f cut = 2/, i.e. all contacting particles whose normal contact force is greater than 
twice the average contact force. Packings for all fi exhibit similar behaviour. 
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D. Contact Network 



The existence of heterogeneous force networks is supported by experimental visualisation. Photoelastic particle 
packings |23| or piles subject to local force perturbations [ p4|j25| , demonstrate inhomogeneity in the magnitude of 
the forces propagating through granular assemblies - "force chains" . However, it is still not clear how relevant these 
structures are in determining the stress state of the system. The 2D simulations of Radjai et al. |ll| suggested 
for compressed granular packings, a distinction can be made between the "strong" force network, those particles in 
contact that carry a force greater than the average normal contact force, and the "weak" force network, the network of 
particles that experience a force smaller than the average. In some theoretical approaches, the strong force network is 
assumed to support all the stress in the system, with the weak force network acting merely as a supporting framework 
to this which can essentially be neglected |2(| . 

To investigate the relative importance of the force networks, we computed the normal force that sub-networks of 
force chains contribute to the bulk average contact normal force. In Fig. |Io|, we varied f cut and then computed the 
fraction of bonds remaining in the force network whose contact force was greater than f cut ('strong' force network), 
and computed the contribution that this network made to the average force. The computation of the relative force 
network contributions in Fig. [l0| indicates only a weak distinction between the 'strong' force network for particle 
contacts with f n > 2, and a weak force network with /„ < 2, say. Therefore it is questionable whether the so-called 
strong network actually does carry most, if not all, of the stress in the system. For example, by going from one curve 
to the other as indicated by the arrow in Fig. 10, wc find 50% of contacts contribute approximately 80% to the average 
contact force. This is a small distinction, and not nearly an order of magnitude difference between the two networks 
that one might expect if the strong forces dominated the weak phase. 




FIG. 10. Frictional packing (fj, = 0.5)-contribution to the bulk average normal contact force and the fraction of particle 
contacts that make up this contribution, as a function of the imposed contact force threshold f C ut- Solid lines are the contribu- 
tions from normal contact forces /„ larger than the threshold f cut and dashed lines are for the forces that are smaller than the 
threshold. Thick solid line: fractional contribution to the average normal contact force for contacts with /„ > f C ut- Thin solid 
line: the fraction of contacts with /„ > f cu t- Thick dashed line: percentage contribution to the force coming from contacts 
with /„ < fcut- Thin dashed line: the fraction of contacts with /„ < f cu t- The arrow indicates the example where 50% of 
particle contacts contribute to 80% of the bulk average contact force. Packings for all /i exhibit similar behaviour. 



A related question is the stability of the relative networks. The concept of fragility [£7j suggests in the limit one can 
clearly distinguish the strong network from the weak network, the strong phase should be a minimally coordinated, 
particle network. For a 3D frictional packing this suggests a coordination number z = 4 p8] |. To calculate the 
network-averaged coordination number of a subset of particles, the contacting neighbours of the chosen network need 
to be included. In Fig. [H] we draw a schematic for determining the coordination number given a sub-network of 
particles (denoted by the grey particles), knowing the list of network neighbours (white particles). 
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FIG. 11. Schematic for computing the coordination number of a subset network of particles. If f cut determines the grey 
particles to belong to the force network, then to compute the coordination number of this network we need to know all contact 
neighbours (grey and white particles) of the given sub-network. 

Computation of the coordination number z for packings with different fi, over a range of cut-off values f cu t is 
shown in Fig. [l2|. The network-averaged coordination number of particle clusters, based on the forces that they carry, 
decreases monotonically from the bulk averaged coordination (f cu t — 0) to approximately z = 1. It appears that 
f cut rj 2, represents some limit in the system in the sense that for f cut > 2, the average size of particle clusters 
contributing are particle pairs, i.e. the largest cluster that propagates large forces is only of size two. 
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FIG. 12. Coordination number z for packings with different /x, of particle networks as a function of the force cut-off f cut that 
determines whether they belong to the network or not. 

Additionally, we provide examples of force network realizations. In Fig. [l^ we show two force network configurations 
of a slice approximately four particle diameters thick taken from the centre of the large, frictional, periodic system 
(N = 128000 and /i = 0.5). We only show bonds whose force is greater than twice the average. Figure [l3](a) is the 
force network for the absolute normal contact forces F without depth normalisation for F cut = 2F. This corresponds 
to a section through the middle of a wide sandpile. Figure |l^(a) can be compared to the 2D experimental realization 
in Refs. |23||29| ]. If we show all bonds, the force network is dense with many weak forces. This may be an indication 
of the relative sensitivity of the experimental visualisation technique which clearly cannot resolve the smallest forces. 
The number of large forces increases with depth giving a clear indication of the propagation of weight down the pile. 

Figure |l^(b) is the force network for the depth-normalized normal contact forces / with f cut = 2f. We find similar 
configurations for the cylindrically confined packing. Therefore this is the equivalent force network for a confined, 
frictional, unloaded packing. Because of the weight of the particles have been renormalized out of the force (mimicking 
walls that support forces), forces of all magnitudes are seen throughout the pack. In both cases, we find that extended 
force-bearing structures exist over a range of length scales, but do not necessarily transmit the largest forces only. 







FIG. 13. Configurations of force networks for (left panel) the absolute forces with F > 2F, and (right) the depth-normalized 
forces with / > 2/. Red-blue continuous colour scheme is used to indicate the relative magnitude of the forces with red 
corresponding to forces closest to the lower threshold and blue are the largest forces. Results for periodic packings with /i = 0.5 
and N = 128000. Networks for different /j, appear very similar. The black frame denotes the size of the simulation cell. 



IV. CONCLUSIONS 

We have demonstrated that large scale simulations of granular packings offer insight into the effects of particle 
friction on measurements of the distribution of particle-particle and particle- wall contact forces P(f). Our detailed 
comparison between simulation and empirical fits, obtained from experiments showed moderate agreement. How- 
ever, we were only able to fit our data over the full range in / after renormalising our data by neglecting the smallest 
forces in the system (using the partial data set). We reason that this is an appropriate way to account for the limited 
resolution in experiment. Wc also reiterate the fact that many simulation and numerical studies of force distributions 
do not show a clear exponential tail either at large / and we believe this may partly be due to the resolution of very 
small forces in such computer experiments that affect the total normalisation parameters. 

We were able to discern the influence that friction plays on P(f) in the small force region. The fraction of particle- 
particle contacts that experience very small forces increases with friction even though the total number of contacts 
decreases with increasing \x ||. Excluded volume effects rather than the functional form of the force law appear to 
dominate the bulk behaviour of the system for dense packings. Our studies of very large systems show that the tails 
of P(f) become marginally steeper with decreasing friction, i.e. j3 in Eq. ^| increases as fi decreases. Our ongoing work 
on simulating confined packings will investigate some of these issues further Jl?j| . 

The force-force spatial distribution function and contact point radial distribution function indicate that spatial 
correlations between the contact forces and the positions of the contact extends out only to approximately three 
particle diameters. This shows that force correlations dissipate quickly in the bulk and that the force transmission 
network propagates locally but becomes diffuse rapidly. On introducing a force cut-off scheme to analyse force-bearing 
structures, we found no clear evidence for distinct 'weak' and 'strong' force phases. 
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In general this discussion was only concerned with unloaded or unperturbed amorphous granular packings. In this 
sense we have provided information on the "reference state" of a granular material from a microstructural point of 
view. This state is rather insensitive to friction and is primarily determined by construction history [^0|. Although 
we have not investigated the perturbed state or response function of these systems |24), it is likely that particle 
properties play a much more significant role in the response of a granular system than in the static state. Some 
theoretical treatments |2?jj3l|] on force chain analysis may benefit from the information of this unperturbed system 
when calculating the resulting response of such a system under some force perturbation. In fact, the contact angle 
distribution in Fig. || seems to suggest that the 'splitting angle' 8 S = 90° — 9, in the language of Re f. Ipj| , does 
seem lie predominantly in the range < 9 S < 30°. Comparing the experimental visualisation in Ref. p4| and the 
theoretical model in Ref. |5l] , the force chain analysis can be thought of as a superposition of force chains on top of 
the background force network. 

This work was supported by the Division of Materials Science and Engineering, Office of Science, U.S. Department 
of Energy. This collaboration was performed under the auspices of the DOE Center of Excellence for the Synthesis and 
Processing of Advanced Materials. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed 
Martin Company, for the United States Department of Energy under Contract DE-AC04-94AL85000. 
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